Towards mechanical entanglement in nano-electromechanical devices 
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We study arrays of mechanical oscillators in the quantum domain and demonstrate how the motions of distant 
oscillators can be entangled without the need for control of individual oscillators and without a direct interaction 
between them. These oscillators are thought of as being members of an array of nano-electromechanical res- 
onators with a voltage being applicable between neighboring resonators. Sudden non-adiabatic switching of the 
interaction results in a squeezing of the states of the mechanical oscillators, leading to an entanglement transport 
in chains of mechanical oscillators. We discuss spatial dimensions, Q-factors, temperatures and decoherence 
sources in some detail, and find a distinct robustness of the entanglement in the canonical coordinates in such a 
scheme. We also briefly discuss the challenging aspect of detection of the generated entanglement. 
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In 1959 Richard Feynman suggested in a famous talk that 
it appears to be a fruitful enterprise to think about manipulat- 
ing and controlling mechanical devices at a very small scale. 
Since then, the study of micro-electromechanical (MEMS) 
and even nano-electromechanical systems (NEMS) has de- 
veloped into a mature field of research 01 lil • Me- 
chanical oscillators with spatial dimensions of a few nanome- 
ters and very high frequencies can now be routinely manufac- 
tured. Applications of such NEMS range from mechanically- 
detected magnetic resonance imaging, sensing of biochemical 
systems, and ultrasensitive probing of thermal transport. In- 
deed, the NEMS devices that are presently manufactured in 
experimental studies are close to or already on the verge of 
the quantum limit While first quantum effects 

are already being observed and studied, it is interesting to see 
to what extent it is feasible to prepare nano-scale mechanical 
oscillators in states where the quantum nature becomes most 
manifest: in states that are genuinely entangled in the canon- 
ical coordinates of position and momentum. This can be in- 
teresting for a variety of reasons. Firstly, it provides another 
stepping stone towards quantum state control and quantum in- 
formation processing in mechanical systems. This is particu- 
larly fascinating as these systems are macroscopic consisting 
of many million atoms. They would therefore also permit the 
exploration of the limiting region between the quantum and 
the classical world. This might be facilitated by another ap- 
plication of entanglement namely its use to enhance quantum 
measurement schemes where entangled states represent a very 
sensitive probe. 

The key question that will be addressed in this letter is how 
it is possible to entangle mechanical oscillators well separated 
in space, without the need for making them interact directly 
and with a minimum need for individual local control which 
is difficult to achieve at the nano-level. This will be accom- 
plished by triggering squeezing and entanglement by a global 
non-adiabatic change of the interaction strength in a linear ar- 
ray of oscillators, but without individually addressing any of 
the oscillators of the array. In this way, one can achieve long- 
range entanglement that will persist over length scales that 



are much larger than the typical entanglement length for the 
ground state of the system |6]. The physics underlying this 
approach, especially the non-adiabaticity requirement, will be 
discussed in more detail lateron. Several schemes to probe 
quantum coherence of mechanical resonators in different se- 
tups and situations have been proposed so far 0,0]. Notably, 
while the earlier proposal of entangling macroscopic oscilla- 
tors |U] entangles two adjacent oscillators in the context of a 
different physical setup, our scheme allows, without the need 
for individual local control, for entanglement in the canoni- 
cal coordinates between non-adjacent (and possibly distant) 
microscopic oscillators by entanglement transport in a chain. 

The setup that we will consider is an array of double- 
clamped coupled nano-mechanical oscillators as has been ex- 
erimentally studied in the micro-mechanical realm in Ref. 
We assume that the beams are arranged in such a man- 
ner that between adjacent oscillators a controlled and tunable 
interaction can be introduced. In Ref. 01 this is experimen- 
tally achieved by applying a voltage between adjacent beams 
made from gold fabricated on a semiconductor membrane that 
are ordered alternatingly. This induces to a good approxima- 
tion a nearest neighbor interaction that can be controlled in 
strength. The oscillators are assumed to be cooled to tempera- 
tures such that kT/huj -C 1 with u> being the fundamental fre- 
quency of the oscillators, such that the array is operated deeply 
in the quantum regime. Before we discuss the time and en- 
ergy scales that would be required to achieve this regime, we 
will exemplify the mechanism, without taking sources of error 
and decoherence mechanisms into account, as we will discuss 
these in some detail later. We start with the Hamiltonian of N 
quantum oscillators of mass m and eigenfrequency u ordered 
on a one-dimensional lattice, with nearest-neighbor interac- 
tion of strength c. Setting h = 1 and using the qu = ql^/mui, 
and pk = p' k /y/muj, where q' k and p' k are the canonical posi- 
tion and momentum of the oscillators we find 

N 

H =^J2{Pk+ + 2c ) - 2 «Mfc+i) . 

fc=l 

For the moment, we assume for simplicity periodic bound- 
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ary conditions, i.e., qN+i = qi, but this requirement will 
be relaxed later, and set u> = 1, as in this ideal treatment 
this merely corresponds to a rescaling of the time scale. The 
normal coordinates are related to the previous ones by a dis- 
crete Fourier transform, q k = (l/\ /r N)^2^ =1 e 2mkl / N Qi, 
p k = (l/\^N)J2fLi e ~ 2 ^ lkl/Np i- In these normal coordi- 
nates, satisfying Q k = Qjy-fc aR d Pk = Pjf-k' tne Hamilto- 
nian can be written in the form 

1 N 

H= oJ2 ( PfcjP fe + ( J + 4csm 2 (TTk/N)y/ 2 Q k Q{) , 

k=l 

annihilation and creation operators, and expressing the time 
dependent operators Qk(t) and Pk(t) in terms of these opera- 
tors, one arrives at the Heisenberg equations of motion for the 
original canonical coordinates 

N 

Qk(t) = 5Z(gr(0)/r-fe(i) +p r (0)g r - k {t)) 
r=l 

andp k (t) = dtqk(t), where we have defined the two functions 

9k{t) = ££Li e 2 ™ kl / N an(fjit)/(Nu t ) and f k (t) = d t g k (t). 
In this paper we are dealing with states that are Gaussian, 
i.e., states whose characteristic function or Wigner function 
is a Gaussian. As such, it is completely characterized by 
the first and second moments frill . The first moments will 
not be directly relevant for our purposes. The second mo- 
ments can be arranged in the symmetric 27V x 27V-covariance 
matrix T R> s = 2Re((i? - {R))(S - (5))), where R and S 
stand for the canonical operators qi,. ■ ■ ,q n and pi, . . . ,p n . 
At this point, we assume that for times t < 0, the oscil- 
lators are not interacting and are in the ground state. This 
implies that T qnqm = T PnPm — 5 n , m , and T qnVm = 0, for 
n,m = 1, . . . , N. 

In the setting of this paper, we will assume for t < the 
interaction is switched off and the system is in its ground state 
and time-independent. At time t = the interaction is then 
switched on instantaneously to ensure non-adiabaticity and 
consequently the system is out of equilibrium and evolving 
in time for t > according to the equations of motion for the 
second moments given by 

r 9„9 m ( < ) = (an,m(t) + d n , m (t))/2, 
r 9 „ Pm (t) = (&n,m(*) + e„,m(i))A 
(t) = (c„ 

,m\J') ~r Q>n,rn (*))/2, 

Where Cl n ,ni = ^2k—l fk—nfk—m: ^n,rn = dtdm.n/^i 
C-n,m — k.—l ^tfk—n^tfk—ma d n ra — ~^2k—l Qk — nQk — m-i 

and e n , m = d t d n , m /2 . 

Before we discuss in detail the non-adiabaticity require- 
ment and other idealizations as well as the physics behind this 
approach we demonstrate the success of the approach. We 
are now in the position to study the entanglement of two very 
distant oscillators when we ignore (trace out) all the others. 
The chain is translationally invariant, and hence, a single os- 
cillator, say labeled 1, can be singled out, and we may look at 
the degree of entanglement as a function of time and discrete 



distance. We quantify the degree of entanglement in terms of 
the logarithmic negativity, defined as En(p) = log \\p TA || i 
for states p, where p TA is the partial transpose and ||.||i de- 
notes the trace-norm. The logarithmic negativity is an upper 
bound for distillable entanglement and has an interpretation 
of an asymptotic preparation cost and bounds the distillable 
entanglement frill . 

Before we consider the entanglement created in this way, 
let us first remind ourselves about the entanglement structure 
of the ground state of the harmonic lattice Hamiltonian: there, 
the bi-partite entanglement between two distinguished oscil- 
lators is only non-zero for nearest neighbors. Next-to-nearest 
neighbors are already separable for all parameters, as are more 
distant oscillators, even in case of an arbitrarily large correla- 
tion length of the chain when approaching criticality, as has 
been demonstrated in Ref . . 
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FIG. 1: The degree of entanglement as a function of time between 
two oscillators in a chain of length 8 with periodic boundary condi- 
tions (left) and open boundary conditions (right). On the right hand 
side, the first and the last oscillator in the chain are considered, on the 
left hand side the two diametrically opposed oscillators. The values 
for c are in the above units c = 0.3, c = 0.2, and c = 0.1 (depicted 
in light, medium, and dark grey). 

This is very much in contrast to the situation encountered 
here: Astonishingly indeed, we find that even very distant os- 
cillators become significantly entangled over time. This de- 
pendence is depicted in Fig.^(for periodic boundary condi- 
tions according to the above formalism, and numerically for 
open boundary conditions). For a time interval [0, to), t.Q > 0, 
the state of the oscillators with labels 1 and n is separable, 
then, for t > to it becomes entangled. This time to is ap- 
proximately given by to ~ n/(27f2), There is what can be 
called a finite 'speed of propagation' of the quantum corre- 
lations, which is in fact closely related to the speed of sound 
in this chain. The amount of entanglement roughly falls off 
as 1/ n, but becomes strictly zero after a finite distance. For 
c = 0.1, for example, this happens for n larger than 500. This 
long-range nature of the entanglement is remarkable indeed. 

The central idea behind the method above is the well-known 
fact that an instantaneous change in the potential of a single 
harmonic oscillator in its ground state will generally make its 
state time dependent and squeezed. In the same way a change 
in the coupling strength between oscillators drives the systems 
away from equilibrium. In the course of the subsequent time 
evolution the squeezing is then transformed into entanglement 
due to the nearest-neighbor coupling. The origin for this is 
the fact that the time evolution is described by a Hamiltonian 
quadratic in the canonical coordinates and therefore has an 
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effect analogous to passive optical elements. It is well-known 
that a beamsplitter which has squeezed states as an input will 
lead to entangled outputs. This process happens continuously 
in the chain. Finally this entanglement propagates, as every 
other excitation, through the chain and can therefore lead to 
entanglement between distant sites. 

In any realistic setting, this switching can not be instanta- 
neous, and an important question is how fast the switching 
process must be in order to generate significant entanglement 
in the canonical coordinates. Fig.[2]depicts the amount of en- 
tanglement in the first maximum when the interaction strength 
is linearly increased over a time interval [0, t']. We find that 
for times t' < 1, any non-zero switching time is unproblem- 
atic (with very similar behavior found in longer chains). This 
is because the change in coupling strength is faster than any 
eigenfrequency in the system, preventing an adiabatic follow- 
ing. For very slow switching, t' S> 1, most entanglement is 
lost because the system can adiabatically follow the parameter 
change and remains approximately in the ground state. 




FIG. 2: Maximum degree of entanglement between the end points 
as a function of the switching time t' in the above units for an open 
chain of length 8 and c = 0.1. The vertical dotted line represents the 
unit frequency of a free oscillator. 

Let us now turn to the discussion of a realisation in NEMS 
of such an array. Presently, NEMS made from SiC have been 
manufactured experimentally with frequencies around 1 — 10 
GHz, with spatial dimensions of the order of 10 nm l^l- flUl . 
Doubly clamped beams have the advantage of higher funda- 
mental frequencies with the same spatial dimensions. The Ce- 
faclors for NEMS of these dimensions achieve values of sig- 
nificantly more thanQ = 10 3 0G3]. Concerning the extent 
to which the ground state can initially be reached, cooling of 
the oscillators to 10 mK seems feasible fl^l using a helium di- 
lution refrigerator Q|, 0] (for the possibility of the equivalent 
of laser cooling to the ground state, see Ref. 1 17]). 

Needless to say, decoherence mechanisms cannot be en- 
tirely avoided in a quantum system so close to macroscopic 
dimensions. After all, Q-factors describe nothing but the cou- 
pling strength to external degrees of freedom beyond our con- 
trol. Most of the dissipation and decoherence is expected to 
be due to the coupling with the degrees of freedom of the 
substrate to which the resonator is connected. Let us now 
specify the decoherence model fl8ll : In the setting described 
here, we are not in the high temperature limit, but close to 
zero temperature. Secondly, we do not have product initial 
conditions: in a realistic setting, the chain and the environ- 
ment are initially not in a completely uncorrected state, but 
rather in the Gibbs state of the coupled joint system, and then 
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FIG. 3: The degree of entanglement under decoherence and for non- 
zero temperature. Shown is again the situation of a chain of length 8 
for periodic (left) and open (right) boundary conditions. In this plot 
(up to rescaling of the time axis, and a quantification of the coupling 
strength c = 0.3, c = 0.2, and c = 0.1 in terms of the fundamental 
frequency), values are chosen that correspond to the Q-factor Q = 
10 3 , fundamental frequency 5 GHz, and temperature of 10 mK. 



driven away from equilibrium fl9ll . We have hence modeled 
the decoherence process by appending local heat baths con- 
sisting of a finite number M of modes to each of the oscil- 
lators with canonical coordinates qj,Pj for k = 1, . . . , M. 
We choose a (discrete) Ohmic spectral density in which case 
the Langevin equation for the Heisenberg picture position be- 
comes the one of classical Brownian motion in the classical 
limit, i.e., the coupling is specified by the interaction Hamil- 



tonian Hj = (u>j(qj ® Sf=i 9j)i wnere &j = j'A/M, where 
A > is a cut-off frequency. This Hamiltonian induces de- 
coherence and dissipation, and the number ( > has in our 
analysis been chosen in such a manner that the energy dissi- 
pation rate reflects exactly the rate 1/Q corresponding to the 
experimentally found Q-factors (see, e.g., Refs. I2l l22|p . With 
this value of (, the initial state before switching on the inter- 
action is then the Gibbs state of the canonical ensemble of the 
whole chain including the appended heat baths. The resulting 
map is nevertheless a Gaussian operation, such that it is suf- 
ficient to know the second moments to specify entanglement 
properties. This model grasps in the simplest possible manner 
the various noise processes ll2lll in NEMS. 

Fig. [3] depicts the behaviour of the degree of entanglement 
for system parameters that are close to those used in actual ex- 
perimental settings. We see that the scheme is surprisingly ro- 
bust against noise processes and non-zero temperatures. Com- 
parably low Q-factors are not particularly harmful given the 
large speed of propagation; yet too high temperatures, turn 
the correlations into merely classical correlations. This effect 
is evidently more harmful for longer chains. Notably, for two 
oscillators, quite large values of the degree of entanglement 
can be achieved. For example, for a two-oscillator system, 
with system parameters as in Fig. [3] the degree of entangle- 
ment as quantified in terms of the log-negativity reaches val- 
ues larger than 0.6 for c = 0.4. Assuming the ability to cool 
to 10 mK, oscillators with fundamental frequencies of 2 GHz 
would be sufficient to generate entanglement. This would be 
the most feasible starting point in such a scheme. 

The most significant technological challenge in an exper- 
imental realization of this scheme (and actually any scheme 
that involves entanglement in the canonical coordinates of 
oscillators at the nanoscale) is the actual detection of en- 



4 



tanglement. We would need to couple the two chosen 
oscillators to canonical coordinate transducers whose out- 
put is proportional to position and momentum, which is 
fed into an amplifier that produces a classical signal l23l . 
What has to be measured with very high sensitivity are 
the second moments of the canonical coordinates q m , q n , 
p m , and p n , i.e., covariance matrix elements. If not all 
entries can be assessed, bounds of the type E^(p) > 
max (0, - log((((<7„ - <7m) 2 )) + ({p n + p m ) 2 ))/2)) may be 
used to estimate the degree of entanglement. If only a posi- 
tion transducer is available, stroboscopic measurements may 
be employed where only two measurements per cycle are per- 
formed (note that p m (t = 0) = q m (t = 7r/4) and posi- 
tion and momentum are interchanging roles with frequency 
oj) l25ll . Alternatively, continuous single-transducer measure- 
ments may be performed which make use of only a position 
transducer and a sinusoidally modulated output 1 25 1 . This 
leaves us with the problem of measuring position and momen- 
tum with great accuracy: conventional optical transducers, as 
they can be employed in MEMS, are not applicable in NEMS, 
but near-field optical sensors or piezoelectric detectors may 
be used Q]. Refs. E3l describe and make use of a bal- 
anced electronic detection scheme of displacement. The most 
promising to date appears to be a capacitive coupling of an 
electrode placed on a resonator to the gate of a single-electron 
transistor, as studied theoretically in Ref. [24] and experimen- 
tally in Refs. 1 3l 12711 . The sensitivity reached in such setups is 
rapidly increasing, and is presently about a factor of 4.3 away 
from the quantum limit of the considered oscillator l27ll . while 
this factor was still about a 100 a year ago, and it is argued 
that with these techniques, the quantum limit could well be 
reached in the not too far future l3l ll6ll24lE7ll . 



Finally, we would like to briefly mention that the chain of 
mechanical oscillators may also be used in principle as a quan- 
tum channel (compare also Ref. O). If one feeds a half of 
a highly entangled two-mode state into a harmonic chain with 
nearest-neighbor interactions, then any oscillator of the chain 
will at some time be entangled with the kept mode. The func- 
tional behavior of the second moments as a function of time 
can be approximated in terms of Bessel functions 12311 . lead- 
ing to a time t\ of the first arrival of entanglement at the n-th 
oscillator of approximately (linear in n) ti « 2n/ (jfle). 

In this letter, we have presented an elementary method of 
entangling mechanical oscillators on the nano-scale which are 
located at macro scopically different locations at the ends of a 
chain, without the need of addressing each of the oscillators in 
the chain. We have introduced the suggested setup formally, 
and have discussed issues of decoherence and measurement. 
As such, the scheme is not yet a fully feasible scheme ready 
for experimental implementation. Yet, it is the hope that this 
letter can point towards significant next steps that could be 
taken when further exploring the quantum domain with nano- 
electromechanical devices. 
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